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Particle-hole continuum in Dirac sea of graphene has an unique window. It has been predicted to support a 
long lived neutral triplet gapless bosoinc mode that disperses over a wide energy range in the whole Brillouin 
zone [Phys. Rev. Lett. 89, 016402 (2002)]. In this work, using a repulsive Hubbard model, we study the fate of 
such collective mode at zero temperature, in a single layer of graphene doped with electrons or holes. Doping 
modifies the particle-hole continuum and creates additional windows. Neutral spin-1 collective mode survives 
and acquires a gap that increases with doping. It is almost dispersion-less inside the new window. We suggest 
that the spin gap signals emergence of short range singlet correlations and a superconducting ground state for 
doped graphene, supporting a recent suggestion of Pathak, et al. I arXiv:0809.0244v l1. 
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I. INTRODUCTION 

Graphene, a single atomic layer of graphite was the first re- 
alization of a two dimensional elemental metallic structur&i-. 
The salient feature in the electronic spectrum of graphene, 
which differentiates it from other materials, is the presence 
of Dirac cones in its dispersion^'"}, which provides a labora- 
tory to test relativistic type phenomena in the ~ eV energy 
scale. The cone-like spectrum protects the system from im- 
purity scattering^, as well as many-body interactions^, both 
in normali^ and superconducting phases^. Because of the 
robust cone like dispersion with a large Fermi velocity, in 
graphene one observes phenomena such as quantum Hall ef- 
fectiS, which appear only at low temperatures and for very 
pure samples in standard 2D electron gas systems, to survive 
even at ambient temperatures. The high mobility of careers in 
graphene at room temperature which is not appreciably differ- 
ent from its value at the liquid-helium temperatur o "V ' is an- 
other promising property of graphene for device applications. 
Stacking the graphene with further layers produces graphene 
multi-layers, such as bilayer, etc. For few layers, due to quan- 
tum size effect, the cone like dispersion is replaced by other 
types of chiral dispersions'*. However, when the number of 
stacks is large enough (~ 10) to approach the bulk limit of 
graphite, the cone like dispersion is again recovered^^, ex- 
cept for small electron-hole pockets at very low-energies ~ 40 
meV. Therefore the phenomena driven by the Dirac nature of 
careers is common in graphene and graphit e 1 3, 14. i s. 16, 17. 1 8. 19 ^ jj. 
should also be shared with the recently fabricated multi-layer 
epitaxial graphene^*^. 

After the pioneering work of Wallace- on the tight binding 
band picture for the electronic structure of pure and undoped 
graphene and graphite, it has become popular to take into ac- 
count the effects of disorder, interactions and doping on top a a 
band picture"*. Nevertheless, there exists an alternate quantum 
chemical approach to the electronic properties of graphene. 
More than half a century ago, Pauling argued that the ground 
state of graphene can be described as a natural extension of 
the resonating valence bond (RVB) state of benzene2ii22i2^, by 
totally ignoring unbound polar (charge fluctuation) configu- 



rations. This overemphasize on neutral configuration makes 
graphene a Mott insulator. But graphene and graphite are both 
semi metals in reality. Since then both band aspect and Paul- 
ing's singlet correlations have been argued to be present in 
graphene and two important consequences have been brought 
theoretically. One of them is existence of gapless neutral 
triplet bosonic mode^^^ for neutral graphene, which was in- 
terpreted as a two-spinon bound excited state in a long-range 
RVB state. The other is a suggestion of high temperature su- 
perconductivity in doped graphena^iS^ and other excotic su- 
perconducting states^^isi. An insulating RVB state is found to 
be stabilizeds^ at least in the Mott insulating side of the phase 
diagram^. There is a lattice gauge theory simulation of 2 + 1 
dimensional QED which predicts the critical value of the "fine 
structure" constant in graphene can be crossed for graphene 
in vacuum2^. In this scenario, the ground state of graphene in 
vacuum is expected to be a Mott insulator. 

After the initial success with fabrication of undoped 
graphene, it has become possible to dope the Dirac sea with 
electrons/holes by applying appropriate gate voltage"*. One 
can also cover a graphene layer with metal atoms in order to 
supply chemical doping of electrons/holes into the Dirac sea 
without significantly modifying their cone-like nature^f . Mo- 
tivated by possibility of doping the Dirac cone, in this work we 
extend our earlier studies to the case of doped Dirac cone. We 
start with the tight-binding model. We employ the RPA ap- 
proximation to analyze the particle-hole scattering processes 
in the triplet channel. The bi-partite nature of honeycomb 
lattice implies that the nearest neighbor tight-binding Hamil- 
tonian considered in our model is particle-hole symmetric. 
Hence doping with electrons and holes are treated on the same 
footing. Therefore here we focus only on the case of electron 
doping. Populating the conduction band of the Dirac cone 
with electrons modifies the particle-hole continuum (PHC) by 
adding a small 2D-like portion and creating small triangular 
windows around the T and K points of the Brillouin zone, 
shown in Fig.[T](c,d,e). This leads to the gap in the dispersion 
of the spin-1 collective mode. The dispersion of the collective 
mode in this window is schematically drawn with red pen in 
Fig- [11(e). The gaped nature of this collective mode in doped 
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graphene in our RVB interpretation translates into the trans- 
formation of long range RVB state of undoped graphens^^ 
to an RVB state with shorter bonds when the Dirac cone is 
doped. 

This paper is organized as follows: In the next section 
we introduce the model and the RPA formulation. Then we 
present the intensity plots of the RPA susceptibility showing 
the dispersion of the neutral triplet collective mode. In the 
small q limit where we use a linearized Dirac theory around 
the K points we obtain closed form expressions for the sus- 
ceptibility. However, the numerical results are valid for the 
whole BZ. We end the paper with a summary and discussions. 



Hence the model we use is the Hubbard model defined as, 



H 



where i,j denote lattice sites of a honeycomb lattice, and a 
stands for the spin of electrons. The spectrum of the band 
limit, [/ = 0, is given by = zte^ corresponding to A = c, w 
for conduction and valence bands, respectively, with. 



y^l + cos(\/3fcy/2)cos(fc:r/2) + 4cos2(fc^/2), (2) 




FIG. 1: (Color online) Particle-hole continuum of graphene band 
structure along important directions of the BriUouin zone, (a) and 
(b) denote the PHC for undoped graphene, while (c) and (d) corre- 
spond to doped graphene. For better comparision, the energies of 
both doped and undoped cases are reported in units of £f. In this 
figure, we have fi — ef = 0.5 eV. Panel (e) is the same as (c) en- 
larged for clarity, which also schematically shows the neutral spin- 1 
collective mode branch. Panel (f) schematically depicts the cone- 
like band picture with two possible sets of inter-band and intra-band 
processes shown by cyan and blue arrows, respectively. As can be 
seen in (f), doping with electrons have two effects: One is to Pauli 
block some of the inter-band transitions which depletes the region 
of momenta around the F and K points, compared to undoped case. 
Second is to add a new 2D like portion to the PHC corresponding 
to intra-band particle-hole excitations denoted by blue in panel (e). 
Hence the triangular windows around F and K points are produced. 



II. FORMULATION OF THE PROBLEM 

As in our earlier work a^^'^^-'^ we study the spin physics 
which involves spin flip (triplet) processes, by focusing on the 
short range part of the interactions on the honeycomb lattice. 



where t ^ 2.8 eV is the hopping amplitude. The lattice pa- 
rameter a is taken as carbon-carbon distance. 

In this model, U is the on-site Coulomb repulsion. Al- 
though the bare value of U in graphene is ~ 4:t — 5t, but within 
the RPA one should not increase U above Uc ^ 2.23^^. The 
underestimation of Uc is a known artifact of RPA, as in the 
sense of Hubbard-Stratonivich transformation, the RPA ap- 
proximation belongs to the family of mean field approxima- 
tions^'. Therefore to be consistent in applying the RPA ap- 
proximation, we restrict ourselves to values of J7 ~ 2t. 

We implement the RPA approximation in the triplet 
particle-hole channel, which is given by^. 



(3) 



Note that the sign of U for triplet and singlet channels is dif- 
ferent. Hence, when the above triplet susceptibility diverges, 
the contribution of the singlet channel to the total susceptibil- 
ity will remain finite. The retarded bare susceptibility x''"'' is 
given by the standard particle-hole form. 



f c,k+q f Xk 



1 



A ^ ^ k-\-q 



(4) 



where N is the number of unit cells, and X = c,v stands for 
conduction and valence bands, respectively. Here ^ is the 
Fermi distribution function, which determines the occupation 
of the state characterized with quantum labels (A, k) and en- 
ergy e^. At T = 0, and for electron doping case correspond- 
ing to /X > 0, the conduction band is partially occupied; i.e. 
from Dirac point to the Fermi level. So there are two types 
of particle-hole excitations: (A) intra-band transition corre- 
sponding to A = c. (B) inter-band transitions corresponding 
to A = w of the above summation. The numerical calcula- 
tion of x^^-* w) for arbitrary qin the BZ is straightforward. 
However, for the low-energy part of the spectrum where the 
Dirac dispersion = hvF\k\ governs the kinetic energy, one 
can obtain closed form formulae for the inter-band and intra- 
band parts of the susceptibility. Therefore in the following 
we further simplify the expressions to proceed with numerical 
evaluation of the integrals, as well as the analytical calcula- 
tions, the details of which is presented in the appendix. 
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A. Intra-band term 

The nature of PHC is shown in Fig. [T] (e). The part of 
PHC shown in blue corresponds to intra-band particle-hole 
processes. As can be seen this part is qualitatively similar to 
the PHC of standard 2D metals with extended Fermi surface. 
The only difference is in the form of dispersion which unlike 
ordinary metals with quadratic dispersion, here one has cone 
like dispersion. For intra-band particle-hole excitations one 
has: 



(0) \ 



In the T ^ where. 



f ^ ^- f ^ 

c^k-\-q c^k 



f - _ 

^ c^k+q 



1 



1 



c,k+q g/3(e£+,--A0 ^ I 

1 



(5) 

(6) 
(7) 



we get. 



(0) ,^ ^ 1 e(/i-£g )-e(/i-eg) 
xinLM = — — — 7-^ — -^^7^- (8) 



^ V k+q k' 



In the first step function we apply a change of variable k+q^ 
—k, and after converting the summation to integral, we obtain: 

(0) X A f ^ e(fcF-fc) 



e{kF - k) 



^H+q-H) 



iO- 



-), 



(9) 



where A = 3\/3a^/2 is the unit cell area. In the analytic 
calculations we use = hvp\k\ which is valid for low ener- 
gies. Therefore in analytic calculations, we will be interested 
in |gj < kp region, where (f is the deviation of momentum 
from r (or K) point. Since the calculation of the above inte- 
gral is standard exercise in mathematical physics, the details 
of the calculation of imaginary and real part of xlntra given 
in the appendix. The final results are given in Equations ( IA6b 
and dASb . These results are valid for small \q\. For arbitrary q 
we evaluate the integrals with numerical quadratures. 



B. Inter-band term 

Inter-band particle-hole excitations in undoped graphene 
can have any energy from zero up to a maximum value deter- 
mined by the shape of band structure. But in the case of doped 
graphene, the inter-band particle-hole excitations of very low 
energy (near F and K points) are Pauli blocked. Hence in 
Fig. [T] the low-energy part of the PHC of undoped graphene 
(Fig.[T]a, b) is depleted to give the PHC of Fig.[T](c,d) which 
is enlarged in (e). For the inter-band processes one has. 



(0) \ 
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N ^ huj 

k 



f - ^~ f 

c,k-\-q ''v,i 



(10) 



which in an analogous way to Eq. (|9]l simplifies to, 

x^lM.-) = ^Jd'kx 

e(fcF - k) 1 



(11) 



^ - i^k + ^k + q) + ^ - (^fe + o + ^k) + »0+ 



k+q k' 



When working with the linearized low-energy theory, the lim- 
its are from to k^, where k^ is the momentum cutoff for 
the linearized theory, defined in such a way that the area of 
a model circular BZ with cone is equal to the original BZ: 
Trfc^ = Abz = 8-\/37r^/(9a^). Substituting the Unear disper- 
sion in polar coordinates, the inter-band part simplifies to. 



Ikp Jo 
1 



-A 



2ir 



hw - hvpik + \/k'^ + g2 + 2kqcos(t)) + iO+ 



(12) 



The imaginary and real part of the above XhJci calculated in 
the appendix. The final results we use in our plots are given by 
Equations (lAlOb and (IA12b which are again valid for values 
of \q\ which are small compared to fc^?. 



III. ANALYTIC AND NUMERICAL RESULTS 

Within the RPA approximation, the dispersion of the neu- 
tral triplet collective mode can be seen as enhancement of the 
susceptibility given in Eq. (|3]l. Since the Xtli^ct ^ retarded 
response satisfying the Kramers-Kronig relation, the imagi- 
nary (dissipative) part is identically zero in regions where the 
denominator of Eq. Q might have zeros: 



0. 



(13) 



As an alternative to solving the above equations, we give the 
contour plots of S^xfrfp^ct- 

To make the comparison of the original theory (in the full 
BZ) with that of the linearized theory (one Dirac cone in a cir- 
cular BZ), one needs to keep in mind that the values of U for 
the renormalized theory should be renormalized to the band- 
width kci'F- Also with the assumption of only one cone in 
the circular BZ, the kinetic energy the particle-hole processes 
occurring at each cone are underestimated by a factor of 2. 
Therefore when we compare a graph from the analytical re- 
sult obtained for the linearized theory with a numerical result, 
one should keep in mind to properly interpret the value of an- 
alytical results according to U ^ 0.4548L/. For clarity, in 
figures related to analytical results, we have performed this 
conversion. 

Fig.|2]shows the intensity plots of RPA triplet susceptibility 
for various values of U, using the analytic results of Equa- 
tions ( |A6] |, dASl l. dAlOb and ( IA12b . The doping is fixed by the 
value of chemical potential /i = 0.4 eV. Fig. [3] on the other 
hand shows the same quantity obtained from numerical eval- 
uation of the x'"' integral. The borders of the window outside 
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FIG. 2: (Color online) Contour plot for jRxtr ipiot u sing the analytic 
formulae obtained in Equations l |A6| l, l |A8l l, l lAlOt and l lA12t . The 
chemical potential is fixed at fi — 0.4 eV. Different panels corre- 
spond to (a) U=1.82t, (b) U=2.00t, (c) U=2.18t. Borders of the trian- 
gular region of Fig.[T](e) are shown in white here. 



FIG. 3: (Color online) Numerical calculation of the spin-1 collective 
mode for (a) U=1.8t, (b) U=2t, (c) U=2.2t and fi = 0.4 eV. Borders 
of the triangular region of Fig.E(e) are shown in white here. 



the particle-hole continuum where the collective mode can ex- 
ist is denoted by white cuts in both figures. These figures 
are qualitatively similar. The precise shapes of the analytical 
and numerical dispersion are not expected to be identical-^. 
Because the numerical calculation takes the details of band 
structure e.g. curvatures and warping into account, while the 
analytical result is a model calculation assuming a circular BZ 
with one cone, where the kinetic energy part is assumed to 
have the form of a cone all over the circular BZ, without any 
curvature or warping. 

In Figs.|4]|5]we plot the dispersion of collective mode using 
analytic and numerical results, respectively. In this figures, 
we change the doping level /i for fixed values of U. Again the 
figures are qualitatively similar. 

As can be seen in Figs.|2]|3]|4]and|5] the neutral triplet col- 
lective mode in the small q region has essentially no disper- 
sion and behaves as Einstein-like mode for the triplet particle- 
hole channel. For smaller values of U in Fig. [3] the collective 
mode is very close to the upper edge and is about to enter 
the PHC. With increasing the value of U, the electron-hole 
pairs are bound tighter into triplet packs, and hence the "bind- 
ing energy" (measured from the upper edge of the PHC) in- 
creases. Near the limiting value of Uc ~ 2.23t, the triplet 
branch still exists in the gapful form and the gap does not van- 
ish. Figures |4] and |5] show that the value of uigap very weakly 
decreases as one increases the fi in the range reported here. 
The behavior at very small /i will be different (see Eq.[T6]l. 

Finally in figure |6] we plot the T — M cut in the BZ of the 
honeycomb lattice. Comparing this result with the one cor- 
responding to undoped case, where the collective mode will 
be gapless at q = Oi^^, in doped case the intra-band levels 



generated in the PHC push the small q part of the dispersion 
higher in energy, thereby making it gapful. The portion added 
to the PHC of undoped system is a 2D Fermi liquid-like part 
inside which the collective mode can not exists and gets Lan- 
dau damped. The Einstein-Uke dispersion in the triangular 
window around the F point exists in the triangular window 
around the K point as well. 

The gaped nature of spin- 1 collective mode in the triangu- 
lar window of the doped graphene is essentially a result of 
repulsion between the blue PHC states (Fig.[T]e). A simple 
asymptotic analysis of the contribution of these 2D like states 
with conic dispersion (unlike the quadratic dispersion of or- 
dinary metals) gives SRxIntra ~ ("^^ ~ h?Vpq'^)^^^^ which is 
valid for small q, and arbitrary oj. When this asymptotic form 
inserted in Eq. ( fTST l clearly gives a gap in the dispersion which 
tends to increase by increasing U. This indicates that the gap 
acquired by the neutral triplet collective mode is a result of 
repulsion of the PHC states arising from intra-band processes. 
However, note that as can be seen in Figs.|2][3] the repulsion 
of the inter-band PHC states (denoted by cyan) dominates and 
the resultant trend of the gap is to decrease with increasing U. 
This should be contrasted to the ID like phase space arising 
in the edge of the undoped PHC region which gives rise to a 
gapless dispersion of the collective mode^. 

Our analytical results valid for the small q description of 
this collective mode enables us to find the behavior of the gap 
in the spectrum of this collective mode in terms the doping 
parameters. To find the dependence of this "spin gap" to pa- 
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FIG. 4: (Color online) Analytic calculation of collective mode for (a) 
fi = 0.2, (b) fi = 0.4, (c) fi = 0.6 eV and U = l.St Borders of the 
triangular region ofFig.H(e) are shown in white here. 
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FIG. 5: (Color online) Numerical calculation of collective mode for 
(a) n = 0.4, (b) n = 0.5, (c) /.i = 0.6 eV and U=2t Borders of the 
triangular region ofFig.E(e) are shown in white here. 



rameters, in q 



limit we have: 

Ah 



(14) 



which gives the following dispersion for the spin-1 collective 
mode: 



^spin 



B 



if. 



4q 
vf 



(15) 



with S = ^ 



hvpkc- As can be seen the 



^ and A = nvpt 
constants are such that the above dispersion can be linearized 
in q, in agreement with the density plots presented earlier. To 
obtain the spin gap magnitude, we set (f = in the above 
formula to obtain: 



21n(^ 



(16) 



Since the doping level ji can be most conveniently controlled 
by gate voltages, this formula can have experimental conse- 
quences. One can see that the above gap vanishes in the /i ^ 
limit corresponding to the undoped case, in agreement with 
our earlier finding'^. For very small values of doping, /^t <C t, 
the characteristic inverse logarithmic dependence of the spin 
gap dominates. This suggests that optical spectroscopies with 
polarized radiation are likely to help in tracing the location of 
bound state in order to examine the peculiar inverse logarith- 
mic ji dependence in our formula. 
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FIG. 6: (Color online) Numerical calculation of collective mode for 
U=2t and /i = 0.4 eV along V — M direction of the Brillouin zone. 
The PHC is more clear in Fig.[T](c). 

IV. SUMMARY AND DISCUSSION 

We investigated the dispersion of the triplet neutral collec- 
tive mode which was predicted to exist in undoped Dirac liq- 
uidsi& in presence of doping. Populating the conduction band 
of the Dirac cone generates new levels in the PHC which re- 
pel the triplet bound state and make it gapful via a genuine in- 
verse logarithmic dependence on the chemical potential. This 
new portion added to the PHC also causes Landau damping of 
some parts of the neutral triplet collective mode as compared 
to the undoped case. 

Such triplet collective modes are characteristic of Dirac liq- 
uids, as opposed to ordinary metallic states with extended 
Fermi surface. In the Dirac cone case, there will be windows 
below the PHC near to F and K points, with no free particle- 
hole pairs allowed. This windows provide the essential phase 
space for the formation of a gaped triplet particle-hole pair. 
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This bound state for the undoped case was interpreted as a 
two-spinon bound state on top of a long range RVB ground 
statai^^. For small values of doping where this RVB pic- 
turo2S can be assumed to prevail by continuity, the gaped na- 
ture of such collective mode can be interpreted as a new RVB 
state with shorter bonds. It would be interesting experimen- 
tal question to search for this inverse-logarithmic dependence 
of the bound state energy. This also suggests a genuine way 
shortening the valence bonds by "screening" the bonds. 

Development of short range singlet correlations in zero mo- 
mentum channel in a conductor is synonymous to develop- 
ment of superconducting correlations. In an independent vari- 
ational study, using the same repulsive Hubbard model that 
we have used, Pathak and collaborators have found supercon- 
ductivity in graphene, for a range of doping. At an optimal 
doping they have predicted high temperature (Kosterlitz type) 
superconductivity. We interpret our spin- 1 collective mode as 
signalling growth of short range spin singlet correlation in the 
ground state. Further, in a self consistent theory, our spin- 
1 collective mode should modify the metallic ground state 
and become the spin- 1 collective mode of the superconducting 
ground state. 



where z = ui/vp- The delta integrals are simpUfied by using 
the formula S{f{x)) ~ J2s \v^'f{x)\ - where s is a root of 
f{x). Let us define, 



f{k,q,(t)) =z±{k- ^/B + +2kqcos(j)). (A2) 
In term of new variable u = cos 0, the root of f{k, q,<j>) = 

is, u = ^-^^j^^, which gives 



V^f{k,q,u)^~-^ (A3) 



Substituting, 
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APPENDIX A: ANALYTIC EVALUATION OF INTEGRALS 
FOR THE LOW-ENERGY THEORY 



du 



(A4) 



in Eq. iAli the u integral becomes trivial and we are left with 
the following integration over the radial variable k: 



Let us calculate the integrals needed in Eqs. (19112b . 



1. Intra-band term 



Using the formula 



+iO+ 



iTr6{x), the imaginary 



part can be most conveniently written as. 



^Xintra(9'2) = 



-A 



dk X 



2'Khvp JQ 

(fc-z)e(fc-^) {k + z)&{k 



(A5) 



q — z > 



2kq 



-q^+2zk \2 
2kq ' 



^Xintral9'^i 



A 

"4^ 



d^k 



27r 



j kdk J d(j)[5{z - k + y/k 
~6{z+ ^ 



^ q^ + 2kq COS < 



\/k^+q^T2k^^cos4))], 



The step functions in integrand correspond to the particle-hole 
(p-h) continuum. We restrict ourselves to the q < 2kF region. 
In this region, the dissipative part of intra-band processes is 
non-zero only when oj < qvp- The radial integration can be 
(A I) performed to give^^i^. 
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^X|°tra('?:^) = 
q 

IGnhvp - z^/q'^ 



1 - 2^ In ' ' 

1 J \\ 1 



2kp — z\^ _^ I ^ — z\ I f 2kp — z 



9 / \ V 9 



For calculation of 3fixintia('?*i ^) ' directly use Eq. ^ and we find: 



2fc_F + z 
9 



2fcF + ^\ f2kF + z 



q 



d^k 



z — k + ^/k"^ + q"^ + 2kq cos z + /s — i/A;^ + f/^ + 2fcg cos ( 



fcdfc / d(l)- 



k — ^k^ + q'^ — 2kq cos < 



2Tr^hvp Jq z'^ - {k- + _ 2kq cos cj))'^ ' 

The (/) integral can be evaluated and simplified to give, 

A I-'"'- 



( k{z~k) 



k{z + k) 



y [(g - fc)2 + {z- fc)2] [((7 + fc)2 _ _ A:)2] ^[(z + fc)2 - (9 - fc)2][(2 + kf - (g + fc)2] 

v4 



1 



(g2_2z2)ln ( 



2fcF + z , 



,2fcF + z. 



1 



.(,2 _ 2z^) In (^^) + AC-hL-±y - 1 - (g^ - 2z2) In ((i) + J(J)2 - 1 



-(g^-2z^)ln(-(-) + J(-)2-lj-g^( 



2 o.2m„ / , //^^2 1^ ^2^2fcF-z^ \^2kp-z^^ ^ ^ ^ ^kp + z ^ l^ 2kp + z ^^ 



(A6) 



(A7) 



-)l (- 



-)i (- 



)2-l 



(A8) 



Demanding the expressions under square root to be positive, gives the following region for the window protected from the 
particle-hole energy levels: q < z < —q + 2kp. 



2. Inter-band term 



Starting from Eq. ( fT2] i. the imaginary part of Xintcr '-^'^ written as 

c- (0) \ A 

First we do integration on cj), and we 



, kdk / 



5 (^z — k — yjk'^ + q^ + 2kq cos 



A 



dk- 



z ~ k 



^nhvpjk, q^l~{ 

A 2z2 - g2 
IQhvp ^z2 _ q2 ' 



2kq 



e(fc 



z + q 



e( 



z-q 



k) 



q + 2kp < z < q + 2kc 



Now we use Kramers-Kronig relation for calculation of ^xlntcr from imaginary part, xl^toT~' 



^xtlM^) 



A 



167rfti;|7(_,+2feF)t,F t^' - w ^^2 -q^vj,' 



(A9) 



(AlO) 



(All) 
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Defining the new variable 77 by the relation tj' = (/Ui? coth(77), the limits of integration 771, 772 are given by coth(772) = \ + 2kc/q, 
coth(7/i) = — 1 + 2A:i?/g, so that we obtain: 



CO (0) N M 



2 ''■02 2coth2(77)-l 
drj— 



A 



-q 



— 2 + 2 



sinh(77) [lu — gw^? coth(7;)] 



2g2(l - 2C.V4) 



2kp 

q 



1-^ 



arctan 



2uj , 



— arctan 



2kF 2ki 



fe_F 



2fe„ 2fc<; 



1+t 



(A12) 



Equations 

(0) , (0) 
A intra ' A inter 



dAlOb and ( IA12b complete the analytic evaluation of the total non-interacting susceptibility 
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